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ral symmetry breaking (DCSB) of QCD vacuum and reproduces meson spectra. 
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5r , produced by phenomenological interaction which consists of a linear confining 

potential, a Coulombic potential, and the 't Hooft determinant interaction. In 
this model, the 't Hooft determinant interaction plays a important role not to 
only 7], 7]' mass difference, but other meson masses through DCSB. 
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1 Introduction 

It is widely accepted that QCD is the correct theory of strong interactions. However, 
it is not so easy to treat the hadron-hadron scattering phenomena by QCD directly. 
Even the lattice QCD may not be able to calculate phase shifts of partial waves in near 
future, although there is an attempt to give some hadron-hadron scattering lengths. 
In this situation it has still meaning to model QCD in a way which has the features of 
QCD to describe hadron-hadron scatterings from the quark freedom. 

According to ref. |1]], we assume the existence of a set of phenomenological inter- 
actions Hpf^^^ such as the residual interaction Hq^j^ — H^f^^^ can be perturbatively 
treated. There exists at the present time no model reproducing the experimental data 
and incorporating simultaneously all the non-perturbative phenomena, such as dynam- 
ical chiral symmetry breaking (DCSB) or confinement, which are the basic aspects of 
QCD. The models which incorporate DCSB and confinement by using the instanta- 
neous confinement potential are at odds with phenomenology [@], |0|, |Q, p, 0. The 
covariant models which preserve the confinement give rather good fitting for the quark 
condensate and the pion decay constant [0], Pj. However, these models must be regu- 
larized for the infrared divergence, and furthermore in the covariant models it is very 
difficult to numerically obtain the quark propagator of the time like region which is 
need to calculate meson masses. 

Because of the difficulty for calculating hadron masses in the covariant way we adopt 
instantaneous interactions which include the confinement potential, the Coulomb po- 
tential with the running coupling constant, and the 't Hooft determinant interaction 
10 as the intermediate range interaction. We show here that the 't Hooft determi- 
nant interaction makes the important contributions to not only t] and rj' meson's mass 
splitting but other meson's masses, because it contributes to DCSB. 

The outline of this paper is the following: In section 2, we find the non-perturbative 
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vacuum using Bogoliubov-Vallatin transformation and write down the Salpeter equa- 
tion which describes meson states. We also comment renormahzation. In section 3, 
we introduce the 't Hooft determinant interaction as a phenomenological interaction 
at intermediate range and show our results. Section4 is devoted to summary. 

2 A phenomenological potential model for vacuum 
and meson spectra 

The Hamiltonian for quarks interacting through an instantaneous fourth component 
potential is given as follows: 



H 



Hq + Hpot 

d^xil)\x){—ioL ■ V + (3m)il){x) 

A" 



+ 2^7"^'^ ''^^y 



i/j'{x) — ip{x) 



V{x - y) 






;i) 



where the row vector ip'^ for three flavour quarks represents [u^ S s^), rh is the current 
mass matrix: 

rh = diag(m„, m^, mj, (2) 

A" 
and — are generators of SU(3) colour group in the fundamental representation. 

The potential is the sum of a linear colour-confinement potential and a Coulombic 

potential with the running coupling constant. In momentum space it is written as 

follows: 

(3) 



V{q) = —r + 



Q 



q 
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where a(g^) is the running coupling constant which varies with momentum transfer q 
According to ref. [0, a(g^) is chosen such as 

r 127r 

I /'QQ O A'l 

I ao (g < go) 



{33-2Nf)\nqyAlcD 



(g > go) 



(4) 



where 

go = AQCZ>exp^^^-|^, (5) 

and Nf is a number of flavour. 

In ordinary perturbation theory, quark fields are expanded in terms of plain wave 
solutions. 

^(^) = E Ij0j^ KP' ^Hp^ ^) + bH-P, -sH-p, -s)] e^P-^, (6) 

where we suppress the colour, flavour and spinor indices. The vacuum state |0) is 
defined as follows: 

aip,s)\0)=b{p,s)\0) = 0. (7) 

Since all other state vectors are constructed on this vacuum, we call it perturbative 
vacuum. 

To find the lower energy state than the perturbative vacuum, we perform the 
Bogoliubov-Vallatin transformation, i.e., quark fields are expanded in terms of chiral 
transformed spinors: 



d^p 



^(^) = E /t^ [«b> ^Mp, s) + P\-P, -s)V{p, s] 



^^P■X 



(27r)3/2 
where 

U{p,s) = exp{j-pep)u{p,s), (9) 

V{p,s) = exp{j-pep)v{-p,-s), (10) 

and 6p is the chiral angle which shows amount of chiral symmetry breaking, a and 
j3 are annihilation operators of pseudo-particle and pseudo-anti-particle, respectively. 
These operators are connected to those of particle as follows: 

a(p, s) = cos 9p ■ a(p, s) — sin 9p { ^ ]-bU—p,—s), (11) 

\ (T ■ p I 

f3\-p,-s) = sin^J'^Q^ ^^_ \.a{p,s) + cos 9p-b\-p,-s). (12) 



Now, we define the non-perturbative vacuum state |0)) , which is the vacuum of 
pseudo-particle and pseudo-anti-particle as 

a{p,s)\0))=f3{p,s)\0)) = 0. (13) 

It is useful to rewrite the Hamiltonian as the normal order with respect to above 
operators. 

H = e+:h: + :Hi: , (14) 

e is the energy of non-perturbative vacuum: 

e = (-6) d^p[msm2{ep + 6p) + p cos 2 (6'p + 5p)] 
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+ 



fd^pd^p'V{p-p') 



(27r)3 

- 7^ Jd^p d^p'Vip - p') sin 2(^p + 5p) sin 2(^p, + V) 

- 7^ jd^pd^p'V{p - p'){p ■ p') cos2(^p + 5p) cos2i9p, + V), (15) 

where the angle 6p which leads to amounts of explicit chiral symmetry breaking because 
of the current mass is defined as 

TTl 

sm26p^^===. (16) 

h is the generalized one body Hamiltonian including self energy effect: 

h = I d^xi)\x){-ioL ■ V + (3m)^{x) 

+ 1 X ]^jd'xdSJ^V{x-y)e'P<''-y^ [ij\x){l - 2^_{p)}iJ{y)\ , (17) 

where A+ and A_ are positive and negative energy projection operators, respectively. 
The angle 6p is determined so that the energy of non-perturbative vacuum may 
be lower than that of perturbative vacuum. Using variation principle, we get the gap 
equation. 

p sin (j)p — m cos (pp 



+ 3 sin 0p 



j -n^y {p-k)p-k cos (pk 
- ^cos(pp j ——^V{p-k)sin(pk 
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= 0, 



where 0p = 2{6p + 5p). 



The procedure stated so far is independent of the particular potential. Since the gap 
equation with the Coulombic potential contains divergent integral, we must renormalize 
wavefunction and mass. The wavefunction and mass renormalization constants are 
given by 



97r Jg^ dq^ [ g^ 

Zm-l = -— / dq q — — , 



(19) 
(20) 



where we must introduce the lower cut-off q^ to avoid the logarithmic divergence. 
Hence the renormalized gap equation is 

p sin 0p + (Z — l)p sin (pp 

— m cos (f)p — {Zm — ^)'m' cos (f)p 

2 f d^k 

+ gSin^pj — — \/(p-fc)p-fccos0fc 



d^k 



- 2 cos vp 
= 0. 



I(2^^^^P-^^'''" 



(21) 



If we solve this renormalized gap equation, we can calculate the quark condensate. 
However, even the quark condensate subtracted the contribution of the explicit chiral 
symmetry breaking is still divergent. To compare with QCD sum rule results for quark 
condensates, it is reasonable to introduce the cut-off Ac=lGeV. In this case the quark 
condensate is calculated as follows, 



(ijij) 



= d^p 

^2^ 



m 



y/p^ + mr 



(22) 



We also make cut-off at f GeV for tfie renormalized gap equation. 

Next step is to construct a meson state on tlie non-perturbative vacuum. We 
consider a meson state wliicli is a bound state of flavour i quark and flavour j anti- 
quark, and postulate that i is not equal to j. Hi is equal to j, we must take the 
wave function of isospin eigenstate. To calculate meson masses, we have to solve the 
Salpeter equation at zero center-of-momentum. It is known that the RPA equation 
is equivalent to the Salpeter equation, but the RPA equations are derived for various 
type of interactions more easily than the Salpeter equation. 

We write down the RPA equations: 

{{0\[H,X]\m)) = {Eo~Ej{{0\X\m)), (23) 

{{0\[H,Y]\m)) = (Eo-ErrMmim)), (24) 

where |0)) and \m)) are the non-perturbative vacuum and the meson state, respectively. 
Inserting the Hamiltonian Eq.(|l]) to above equations, we obtain the following equations, 

[E,{k) + E,{k)]X{k) - ^jj^V{k - k')[X{k') + Y{k')] = MX{k), (25) 



4 r d'^k' 
mk) + E,{k)]Y{k) --j T^,V{k - k')[X{k') + Y{k')] = -MY{k), (26) 



d^k' 
where the forward and backward amplitudes of a meson state are defined as 



X{k) = {{0\X\m)) = Y: U\-k) ((0|/5;(fc)<(-fc)|m)) V'^\-k), (27) 

Y{k) = ((0|F|m))= Y. V\-k){{0\af{-k)f3fik)\m))W\-k), (28) 

M is the mass of the meson state, and Ei and Ej are single particle energy of the 
particle and the anti-particle, respectively. In general above equations depend on the 
center of mass momentum P and the relative momentum k. Here we set P = 0. Dirac 
spinor U{k) and V{k) are eigen-vector of the generalized one body Hamiltonian h in 



momentum space: 

h{k)U{k) = E{k)U{k), (29) 

h{k)V{k) = -E{k)V{k). (30) 

The forward and backward amplitudes are connected to the Salpeter amphtude x(^) 
as follows: 

X{k) = AV(fc)x(fc)Ai(fc) + Ai(fc)x(fc)Ai(fc) 

= X{k) + Y{k), (31) 

Hence the RPA equations are converted to the Salpeter equation: 

[E,{k) + E^{k)] [A\{k)x{k) - x{k)A\{k) 

4 r d^k' r 1 

" 3 1^^^^ ~ ^'^ K(fc)x(fcO - x(fc')A^+(fc)] = Mx(fc), (32) 

and the Salpeter amplitude x obeys the constraint: 

AV(fc)x(fc)-x(fc)Ai(fc) = 0. (33) 

To reduce the Salpeter equation, we expand the Salpeter amplitude in terms of a 
complete set of Dirac matrices [§ 

X{k) = Co{k) + ^ A(fc)p, + A/'o(fc) .a + Y.M{k)- p^a. (34) 

The constraint Eq.(p3D demands 

C^{k) = i^-^^M^{k)-k, (35) 

COS0_ 

C,{k) = -^^M(fc).fc, (36) 

At,(k) = -^^C,(k)k-t^^J^,{k)xk, (37) 

sm 0+ cos (p- 

A7-o(fc) = t^^^C2{k)k-^^^At2ik)xk. (38) 

cos 0_ sm 0+ 



The general form of x satisfying the constraint becomes 



+ M- 
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cos/ 
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. sm 



smi 
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COS( 



cos </>+ / ^ .- . . sm (/)_ " 

P20- + — (cr X fc) + I —k 

sm c)+ cos ©_ 



The part of A^x' ^ x'A.+ in the left-hand side of Eq.(|32D is written as follows: 



cos <p+ ^ - 
pi ^-—p3cr ■ k 



p2 + i- 



smc/; 
.sine 

COS( 



AX{k)xik') - x{k')Ai{k) 
= l—i sin </)_|_ cos (P-C2 — sin 0+ sin 0_ (A/q ■ fc) [ 

+ I i sin 0_|_ cos 0_£i — i cos 0_|_ cos 0_ (A/3 ■ fc) !• 

+ I — i sin 0+ cos 0-A/2 — sin 0+ sin c^^Cji — i cos 0+ cos 0_ (A/q x fc) i 

COS0 



-a • fc 



picr 



smi 



-- . sm (/)_ ^ - 

-psfe + I —{P30- X fe) 

cos 0_ 



+ 1 2 sin 0+ cos 0-A/'i — z cos 0+ cos (p-C'-^k + sin 0+ sin 0_ (A/'3 x k)> 



P2<? + 



cos I 



smi 



, ^ ' , , sm 0_ j- 
(a X k) +t —k 



cos I 
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')i ± 0j 
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and the functions with the prime depend on k'. 

Then the Salpeter equation is reduced to coupled equations, 

Ei + Ej sin I 



-A 



COSi 



+ - lv7r^V{k-k') 
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sin 0, cos 6M + sin 0, sin 0_ ^^Z:'. fc ■ fc 



cos- 



;,/ '^2 



COS 



i sin 0+ sin 0_ ^ — -j- (A/'2 x k ) ■ k 



sm 



M 
2i 



/:i, 



(39) 



(40) 



(41) 



(42) 



Ei + Ej cos (/)_ 

2 sin (/)+ 

2 /• d^k' 
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(2vr) 



y(fc - k') 



sin04. cos(/)_£']^ 



COS( 



+ i cos 0+ cos I 



Sllli 



cos 



-^{Atlxk)-k 



cos 04_ COS ( 



2i ' 



SlIli 



~ JL^i r\i ' K 



(43) 



Ei + £^-,- 



coscp- 
sin(/)_i_ 



■M 



cos 0+ COS ( 



sm 0+ sm (p- 
cosd)- 



{M2 ■ k)k 



2 r d-'k 



3lJ 



3 7 (27r) 



sm 01 sm ( 



:V{k-k') 



— * cos —* ^ I ^ 

sin 0_i_ cos <^-M'<2, + cos 0+ cos 0_ -^ (A/'n X k ) X k 

sm (I) , 



""'^-(A/-,'.fc)fc-.cos0,cos^ "" 



/v' 



JL^tylx, X K 



COS( 



cos- 



;,/ '^2 



M 
2i 



M, (44) 



Ei + -Ej 



d^fc' 



sm( 



COS( 



-M 



cos 04- COS I 



smi 



sm 0+ sm 



cos I 



(M ■ fc)fc 



37 (27r)3 



y(fc-fc') 



cos I 



sin 0+ cos ^-A/"! + cos 0+ cos 0_ -p {N'[ ■ k)k 

sm0L, 



sm 0_i_ sm ( 



sin 0'_ , ,-.., i- ^ ,- , , cos 0', ^, - / 
— (A/^ X fe ) X fc + zsm0+sm0_^ T'^i''' ^^ 



COS( 



smi 



M 
27' 



■Mo. 
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If we make i equal to j, we obtain the same equations as those derived for one flavour 
quark by Yaouanc et, al. [Q] 

The pion decay constant /tt is gotten by using the pion wave function for the center 
of mass momentum P 7^ of pion according to the method of ref . [§ as follows: 

r-Ac d^p 1 2Ci ^^''^ 



U 



(27r)3M^sin, 



(46) 



The Gell-Mann-Oakes-Renner relation |Tl[ is given as 



-2Z^m^{^ij)o = M^f^, (47) 

where {4''4')o is the condensate which is calculated by setting niu = 0. 

Let us discuss numerical results. Quark masses, strength of linear potential, and 
Aqcd are fixed as m^ = m^ = 5MeV, m^ = IQOMeV, a = 4.69fm~^, and Aqcd = 
200MeV/c. The strength of the Coulombic potential ao and the lower cut-off mo- 
mentum ql included in the equations for Z — 1 and Zm — 1 are freely varied so as 
to fit to empirical values of the pion decay constant, and condensation of u and s 
quarks. For oq = 5.05 and ql = 0.86go we get /^ = 93MeV, (uu) = -(263MeV)3 
and (ss) = — (240MeV)^. Factors of counter terms, Z — \ and Zm — 1 are —0.580 
and —0.759, respectively. Results of calculated meson masses are shown in Table 2. 
The pion mass is about one half of the experimental value. Also, other meson masses 
are about 200MeV lower than experimental ones. Single particle energies of quark 
and anti-quark are very large because of the infrared properties of the linear potential. 
Actually, those are divergent at the momentum p = 0. However, large parts of these 
single particle energies are canceled out by the interaction between the quark and the 
anti-quark. In the present model the effect of cancelation is too large. Also, the abso- 
lute value of factor Zm — 1 is too large, so pion and kaon masses become much smaller. 
The reason why the absolute value of factor Zm — 1 is large is that we need to choose 
the large coupling constant ao of the Coulombic interaction. 

It is necessary to introduce such an intermediate range interaction as does not 
essentially contribute the quark-anti-quark interaction part for vector mesons although 
it increases the dynamical mass of single quark. The absolute value of factor Zm — 1 
reduces because the coupling constant Oq of the Coulombic interaction can be chosen 
to be small by introducing this interaction. A candidate of this kind of interaction is 
the 't Hooft determinant interaction which mixes vacua of different winding number. 
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So, let us add this interaction to the Hamiltonian. 



3 The 't Hooft determinant interaction 



As a phenomenological interaction of intermediate range, we introduce the 't Hooft 
determinant interaction: 

Haet = K Jd^x {det [^i{x){l - i5)^j{x)] + det [V^,(a;)(l + 75)^,(2;)] } , (48) 

where the determinant is over flavour indices. 

We apply the same procedure as previous section to the Hamiltonian with the 
't Hooft determinant interaction: 



H = Hq + Hpot + Hdet- 



(49) 



The vacuum energy is 



/ 



where e/ is flavour independent term and 

Since the 't Hooft determinant interaction mixes quark flavours, the gap equation of 
6p and 6** becomes coupled equation: 

p sin (j)p — m cos 0p 
+ gSin^py— -— \/(p-fc)p-fccos0fc 



— - cos 
3 



^'li^^^p-^^''''' 



d^k . f d^k 



/ft n r (i k' 

— -sin0fey— -sin 



(27 



(52) 



p sm (pp — m cos ( 
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2 /■ d^k 

+ ^sm(j)p J -——V{p- k)p-k cos (j)l 

- -cos(plj—-^V{p-k)sin(Pl 



- 8N,{2N^ + 3N, + 1) cos 0^ I J 



(2^ 



■ sm (pk 



0, 



(53) 



where 



^, ^ 2(^« + <5,-), 
p = 2(^; + <5;), 



(54) 
(55) 



and we set m^ = m^^, then the gap equation of 9^ is the same equation as that of 9p. 
We adopt isospin eigenstates for meson states as follows: 

1 



X 



.1=1 



,ud 



X 



1=0 



X 



V2 



V2 



(x^^'-x'') , x"^ 



(x"" + x'') , X' 



(56) 
(57) 



The Salpeter equation for isospin 1=1 and states are given by 



H{k),x\k) 



4 r d^k' 



3 J (27r)J 



iiii\ 



V{k-k') A+{k),x\k 



+ 



+ 



''+'2{N, + 1)K l£^sm 



(27r)3 ^PJ (27r)3 



d^k' 



^[A4fc),7°ra'(fc')7°r, 



i=± 



\^+i 



2N,{N, + 1)K 



d^p f d^k' 

- — — sin <!)„ / - — — 
(27r)3 ^Pj (27r)3 



^Tr(x^(fc')7°r.)[A+(fc),7°r 



i=± 



5ifl 2V2{N, + l)ir 



d^p 



■smi 



rf^A;' 



- (5/,o 2V2Nc{Nc + l)K I t:-^ sin 



(2vr)3— 7(2vr)3;^^ 



(27r)3""""^7 (27r^3 
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5:Tr(x^^(fc')7°r.)[A+(fc),7°r, 



i=± 



Mx'(fc), 



(5^ 



H^ik),x'^{k) 



4 /• rf^A;' 



3 J (27r) 
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Vik-k') AUk),x'%k 
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- 5,,o2v^iV,(iV, + l)i^|^sin0,|^^Tr(x^(fc')7°r.) [A;(fc),7°r, 
= Mx'%k). (59) 

The Salpeter equation for isospin / = 1/2 state is also given by: 

(Ei(fc) + E,{k)) {A^^\k)x{k) - x{k)A^'\k)) 
- I J0^sy(k - k') [A?(fc)x(fc') - x{k')At\k)_ 

+ 2{N, + i)KJ^^ sin0f /|^ E [A«(fc)/r.x(fc')/r. - /r.x(fc')/r.A?(fc) 

+ 2iV,(iV,+ l)ir/^sin0f /^5:Tr(x(fc')/r.) [A«(fc)7°r, - 7°r.Af (fc/ 
= Mx(fe), (60) 

where T± = 1 ± 75. 

The gap equation and the Salpeter equation including 't Hooft determinant inter- 

2/2 

action have to be solved. We introduce a kind of cut-off factor e~^ '^0 to the momentum 
integral parts in which the 't Hooft determinant interaction appears. We choose the 
value of parameter po around 0.6GeV/c because the 't Hooft interaction is adopted 
as the intermediate interaction. We will show results of two cases of po = 0.6GeV/c 
and 0.8GeV/c. In each case, free parameters are the coupling constant of the 't Hooft 
determinant interaction K, the coupling constant of the Coulombic interaction ao, and 
the cut-off Pl. Those are determined so as to fit the pion decay constant and quark 
condensates. 

Including the intermediate interaction, the dynamical mass of quark increases. 
Hence all meson masses increase and the mass of kaon , especially, increases ISOMeV. 
However the calculated pion mass is underestimated 30MeV to the experimental value, 
since mass renormalization is more effective than overestimate of quark condensate. 
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The results are summarized in table 1 and table 2. For meson masses, the difference 
between calculated values and experimental values are within 70MeV. 

4 Summary 

In this paper, we constructed a phenomenological model which reproduces the meson 
spectra. At first, we considered the phenomenological potential model which includes 
the long-range linear potential and the short-range Coulombic potential. The meson 
masses were underestimated because the effect of mass renormalization was too large 
due to the large coupling constant of the Coulomb interaction which is needed to 
reproduce u and s quark condensate and the pion decay constant, and large parts of 
dynamical quark masses were canceled out by the residual interaction energy. Hence 
we introduced the 't Hooft determinant interaction as the intermediate interaction to 
increase the effective dynamical mass of quarks without contributing to the residual 
interaction of vector mesons. 

The 't Hooft determinant interaction played a important role not to only rj, rj' mass 
difference, but other meson masses through DCSB. Setting parameters so as to fit quark 
condensates and the pion decay constant, the difference between calculated masses and 
experimental masses could be within 70MeV. 

Adding the gluon current to the quark current in the Hamiltonian, according to 
ref. [^, we can calculate the glueball spectra in the same method. 

To describe hadron-hadron scattering processes in this model, we need to describe 
moving meson. However it is too difficult to solve the Salpeter equation for non-zero 
center of mass momentum without any approximation. Hence we must make an effort 
to find an appropriate approximation. 
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Table 1 







(uu) [MeV^ 


(ss) [MeV3; 


U [MeV] 


Z^-1 


Z-\ 


Empirical 




-(230 ±25)3 


-(200 ±30)3 


93 






ao = 5.05 

qi = 0.86go 

K = 




-(263)3 


-(240)3 


93 


-0.759 


-0.580 


ao = 2.3 
qi = O.Sqo 
K = 0.0485, po = 


--4.0 


-(283)3 


-(219)3 


92 


-0.605 


-0.344 


ao = 2.4 
qL = 0.84go 
K = 0.173, po = 


3.0 


-(282)3 


-(227)3 


92 


-0.568 


-0.329 



Table 2 







TT 


P 


K 


K* 


V 


r^' 


UJ 





Exp. 




140 


770 


498 


892 


547 


958 


782 


1020 


tto = 5.05 

qL = 0.86go 

K = 




71 


557 


271 


716 








823 


ao = 2.3 
gi = 0.8go 
K = 0.0485, po = 


= 4.0 


103 


769 


420 


867 


473 


930 


774 


951 


ao = 2.4 
qL = 0.84go 
K = 0.173, po = 


3.0 


107 


790 


438 


915 


494 


1006 


796 


1017 
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